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Abstract. We prove the SBV regularity of the characteristic speed of the scalar hyperbolic conservation 
law and SBV-like regularity of the eigenvalue functions of the Jacobian matrix of flux function for general 
systems of conservation laws. 

More precisely, for the equation 

ut + f{u)x=0, « : R+ X R -i> Q C M^, 

we only assume the flux / is function in the scalar case {N = 1) and Jacobian matrix Df has distinct 
real eigenvalues in the system case (N > 2). Using the modification of the main decay estimate in [7] 
and localization method applied in [13], we show that for the scalar equation f'{u) belongs to SBV, and 
for system of conservation laws the scalar measure 

{Du\i{u) ■ ri{u)) {li{u) ■ Ux) 

has no Cantor part, where A^, r^, li are the i-th eigenvalue, i-th right eigenvector and i-th left eigenvector 
of the matrix Df. 
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1. Introduction 

The study of the regularity of solutions to a general system of hyperbolic system of conservation laws 
(1.1) ut + /(u):,=0, u:M+xM^ficM^ 

with initial data 



(1.2) 



u(t = 0) = uq G 



,0) 



is an important topic in the study of hyperbolic equations. In particular, recently there have been 
interesting advances in the analysis of the structure of the measure derivative Dxu{t) of a BV solution to 
genuinely nonlinear scalar equations and hyperbolic systems. The results obtained are that, in addition 
to the BV bounds, the solution enjoys the strong regularity property that no Cantor part in the space 
derivative of u{t) appears out of a countable set of times [1, 7, 13]: the fact that the measure Dxu{t) has 
only absolutely continuous and jump part yields by definition that u{t) £ SBV. 

The main idea of the proof is to find a bounded functional, which is monotonically decreasing in time: 
then one shows that at each time a Cantor part appears the functional has a jump downward, and hence 
one concludes that the SBV regularity of u outside a countable set of times. 

This paper concerns the extension of the results of [7] to the case where the system is only strictly 
hyperbolic, i.e. no assumption on the nonlinear structure of the eigenvalues Ai of Df is done. Clearly, by 
just considering a linearly degenerate eigenvalue, it is fairly easy to sec that the solution u itself cannot 
be in SBV, so the regularity concerns some nonlinear function of u. 

We state the main theorem of this paper: in the following a BV function on R will be considered 
defined everywhere by taking the right continuous derivative. 

Theorem 1.1. Let u be a vanishing viscosity solution of the Cauchy problem for the strictly hyperbolic 
system (3.1) with small BV norm. Then there exists an at most countable set S C M"*" such that the 
measure 

{Du\i{u) ■ ri{u)) (li{u) ■ u^) 
has no Cantor part for every t G M.'^ \ S and i G {1, 2, . . . , N}. 

In the scalar case the above theorem can be rewritten as 
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t : scalar 



Theorem 1.2. Suppose thatu £ BV(M+ x M) is an entropy solution of the scalar conservation law (2.1). 
Then there exists a countable set S C such that for every t e \ S the following holds: 

f{u{t)) G SBVioc(K). 

Since in the genuinely nonlinear case m H- {u) is invertible along the i-th admissible curves [u] (see 
Theorem 3.2 for the definition), it follows that Theorem 4.1 is an extension of the results contained in 
[7] (and Theorem 1.2 is an extension of the results contained in [13] when the source is 0). The example 
contained in Remark 7.2 shows that the results are sharp. 

The main point of the paper is the fact that the wave-tracking approximation for the waves of a 
genuinely nonlinear family does not essentially differ from the wavefront approximations of genuinely 
nonlinear systems: in other words, the wave pattern of a genuinely nonlinear characteristic family for 
a (approximate) solution in a general hyperbolic system has the same structure as if all characteristic 
families are genuinely nonlinear. Thus the analysis carried out in [7] holds also in this case. 

The proof of the above two theorems is done as follows. To introduce the argument in the easiest 
setting, in Section 2, we give a proof for the SBV regularity of the characteristic speed for the general 
scalar conservation laws. The proof is just a slight modification of the proof of Theorem 1.1 in [13]. 

As one sees in the proof of Theorem 1.2, the main tool is to obtain the SBV regularity when only 
one characteristic field is genuinely nonlinear (Corollary 4.2). By inspection, the analysis of [7] relics on 
the wave- front tracking approximation of [8], which assumes that all characteristic fields are genuinely 
nonlinear or linearly degenerate. Thus we devote Sections 3, Section 5.1 to introduce the wave-front 
tracking approximation for general systems [3]. 

The focus of Section 5.2 is the observation that the convergence and regularity estimates of Theorem 
10.4 of [8] still holds for the z-th component of Ux, under the only assumption that the i-th characteristic 
field is genuinely nonlinear: these estimates are needed in order to define the i-th (ei, eo)-shocks and to 
pass to the limit the estimates concerning the interaction, cancellation and jump measures. The latter is 
responsible for the functional controlling the SBV regularity. Theorem 4.1. 

After these estimates, for completeness we repeat the proof of the decay of negative waves in Section 
6.2. Finally we show how to adapt the strategy of the scalar case in Section 7. 
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2. The scalar case 

In this section, we restrict our attention to the scalar conservation laws and motivate our general 
strategy with this comparatively simpler situation. Let us consider the entropy solution to the hyperbolic 
conservation law in one space dimension 



(2.1) 



ut + f{u)x=Q u : M+ X M ^ C M, / e C2(rj, 
U|t=o = Wo uoeBV(M, il). 



It is easy to generalize the SBV regularity result from the convex fiux case to the concave case in the 
following sense. 

Lemma 2.1. [13] Suppose f G C^(M) and \f"{u)\ > 0. Let u G L°°(]R) he an entropy solution of the 
scalar conservation law (2.1). Then exists a countable set S C R such that for every t G M+ \ S the 
following holds: 

u(t,-) G SBVioc(M). 

Further, by Volpert's Chain Rule (Theorem 3.99 of [2]), it follows that /'(u(r, •)) G SBVioc(IK) for 
T G M+ \ S: actually, since /" 7^ 0, the two conditions f'{u{T)) G SBVioc and u{t) G SBVioc are 
equivalent. 

Following the same argument together with the analysis in [13], we can get a SBV regularity of the 
slope of characteristics for the scalar conservation law with general flux. 



Proof of Theorem 1.2. Recall that if w G BV(R+ x M) is an entropy solution, then by the theory of 
entropy solutions it follows that Ut(-) := "(''', ■) G BV(R) for all t G K+. 
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Define the sets 

Jr := {x eM.: u{t,x—) ^ u{t,x+)^, 

Fr:={xeR: /"(it(T,x)) = O}, 

C := {(t,$) e M+ X M : ^ G UF^}. 

Set also Cr '■= Jr U Fr as the r-section of C. 

Since the Cantor part D'^Ur of Du^ and tlie jump part D^'^Ur of Dur are mutually singular, then 
iD'^UrliJr) = 0. Using the fact that ,f"{ur) = on Fr, by Volpert's Chain Rule one obtains 

\D'^f{Ur)\{Cr) < \D'fiUr)\iJr) + \D'=f(Ur)\(Fr) 

= \f"{Ur)D'Ur\{Jr) + \f'{Ur)D'Ur\{Fr) = 0. 

Let {to, xo) G X M\C. Using the finite speed of propagation and the maximum principle for entropy 
solutions and the fact that utg is continuous at xq by the definition of C, it is possible to find a triangle 
of the form 



T{to,xo) ■.= {{t,x) : \x-xo\ <bo-X{t-to), < t - io < 6o/a} 



such that f"{u{t,x)) > cq > 0, for any (t,x) G T{tQ,XQ). Here cq depends on {tQ,xo) and A is the 
maximal speed of propagation, which depends only on the L°°-bound of ut^ (and hence only depends on 
the L°°-bound of u by maximal principle). 

In particular, in T(to, xq) the solution u of (2.1) coincides with the solution of the following problem 



Wt 

w{to,x 



fH, = 0, 
J "to (a;) 



\x -xo\ < bo, 
2^ Ixo-b° "to iy)dy \x-xo\>bo. 



By Lemma 2.1, w{t,-) is SBV regular for any t > out of a countable set of times S{to,xo). Write 
Tr{to,Xo) ■■= T(to,xo) n {i r}, thus Ur^T^(to,xo} and f {ur)^T^(to,xo) are SBV for r €.]to, to + 
b/~X[\S{to,xo). 

Let B be the set of all points of M.'^ x M.\C which are contained in at least one of these triangles. 
(Notice that T{to, xo) is a open set and does not contain the point {to, xo)-) Let C" := M"*" x K \ (B U C). 
We claim that the set Sc' {t G : {t = t} n C" 7^ 0} is at most countable. Indeed, it is enough to 
prove that the set Sk '■= {t G : {t — t} D C O K ^ (d} is a,t most countable for every compact set 
K C M"*" X R when the triangles T{t',x') have a base of fixed length for every {f ,x') G C": it is fairly 
simple to see that in this case the set Sk is finite since (t', x') can not be contained in any other T{t" , x") 
iort' ^t" and{t",x") eC 

Finally, let {T(ti, Xi)}igN be a countable subfamily of the triangles covering B. From the previous 
observation on the function u\_x(ti,xi): the set 

Si := {t : UT^T^{u.x,)t SBY{Tr{ti,Xi))] 
is at most countable. For any t not in the countable set 

Sc U U 5„ 

one obtains the following inequality: 

(2.3) \D-f'{Ur){m < \D-f'{Ur)\({jTr{t„X,)\+\D-f'{Ur)\{Cr) = Q. 

This concludes the proof. □ 
By a standard argument in the theory of BV functions, we have the following result. 



X M) be an entropy solution of the scalar conservation law (2.1). Then 
The difference is that now the function J'{u) is considered as a function of two variable. 



Corollary 2.2. Let ueL°°{ 
f'{u) G SBVioc(K+ X M). 
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Proof. The starting point is that up to a countable set of times, Df'{u{t, •)) has no Cantor part (Theorem 
1.2). From the slicing theory of BV function (Theorem 3.107-108 of [2]), we know that the Cantor part 
of the 2-dimensional measure Dxf'{u) is the integral with respect of t of the Cantor part of Df'{u{t, •)). 
This concludes that Dxf'{u) has no Cantor part, i.e. D%f'{u) = 0. 

By combining Volpert's Chain Rule and the conservation law (2.1), one has 

Dlu = ~nu)Dlu. 

Using Volpert's rule once again, one obtains 

D'inu) = nu)Diu = f"{u)r{u)Diu = nu)D%nu) = o, 

which concludes that also Dtf{u) has no Cantor part. □ 



s:gpn 



3. Notations and settings for general systems 

Throughout the rest of the paper, the symbol 0(1) always denotes a quantity uniformly bounded by 
a constant depending only on the system (3.1). 



e :basic3 



lambda 



assumponri 



■eimon orient 



3.1. Preliminary notation. Consider the Cauchy problem 

iut + f{u)^ = u : M+ X M ^ C K^, / e C2(S1, R), 
\u\t=o = uo woeBV(M,ri). 



(3.1) 



The only assumption is strict hyperbolicity in il; the eigenvalues {Ai(it)}^j^ of the Jacobi matrix 
A{u) = Df{u) satisfy 

Ai(u) < • ■ • < u G fl. 

Furthermore, as we only consider the solutions with small variation, it is not restrictive to assume 
compact. Hence there exist constants {^j}jLo, such that 



(3.2) 



Afc_i < Afe(it) < Afe, yuen, k = l,...,N. 



Let {ri(u)}lLi and {'j(u)}wLi be a base of right and left eigenvectors, depending smoothly on u, such 



that 
(3.3) 



lj{u) ■ ri{u) = 5ij and \n{u)\ = 1, i = I,. . .,N. 



Definition 3.1. For i — l,...,iV, we say that the i-th characteristic field (or i-th family) is genuinely 
nonlinear if 

VAi(u) • rj(w) 7^ for all u e fl, 

and we say that the i-th characteristic field (or i-th family) is linearly degenerate if instead 

VA.i(M) • ri{u) = for all u e fl. 

In the following, if the i-th characteristic field is genuinely nonlinear, instead of (3.3) we normalize 
ri {u) such that 



(3.4) 



VAi(w) • ri{u) = 1. 



In [6], it is proved that if the total variation of uq is sufficiently small, the solutions of the viscous 
parabolic approximation equations 

jut + f{u)r, = euxx, 
|u(0,x) = uo{x), 

arc uniformly bounded, and the limit of as e is called vanishing viscosity solution of (3.1) and it 
is BV function. 
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3.2. Construction of solutions to Riemann problem. The Riemann problem is the Cauchy problem 
(3.1) with piecewise constant initial data of the form 



(3.5) 



Uq 



,R 



X <0, 

X > 0. 



The solution to this problem is the key ingredient for building the front-tracking approximate solution: 
the basic step is the construction of the admissible elementary curve of the k-th family for any give left 
state u^. 

A working definition of admissible elementary curves can be given by means of the following theorem. 



,N, satisfying lims^a l^^sM 'Tkiu), 
T < s and k = 1,...,N, satisfying 



Theorem 3.2. [5, 6] For every u £ ^1, there exist 

(1) N Lipschitz continuous curves s t-^ T!^[u] £ fl, k ~ 1, . 

(2) and N Lipschitz functions (s, r) H- (T^[u](t), with < 
T t— >■ a'^[u\{T) increasing and (Tq[u](0) = Xk{u), 

with the following properties. 

When £ 17, = T^[u^], for some s sufficiently small, the unique vanishing viscosity solution of the 
Riemann problem (3.1)-(3.5) is defined a.e. by 

^ x/t<a^[u^m, 
u{t,x) ■= { T'^[u^] x/t = cr^-[wL](r),r £ [0,s], 
^' x/t> a^^[v^]{s). 

Remark 3.3. If i-th family is genuinely nonlinear, then the Lipschitz curve T"* [t/] Ccin be written cis 



TIM 



R,[u]{s) 
SAu]{s) 



s > 0, 
s < 0, 



where Si[u\ are respectively the rarefaction curve and the Rankine-Hugoniot curve of the i-th 

family with any given point m in J7. And certain elementary weak solution, called rarefaction waves and 
shock waves can be defined along the rarefaction curve and Rankine-Hugoniot curve, for example see [8]. 
The elementary curve Tl[u] is parametrized by 



(3.6) 



h{u)-{Tl[u]-u) 



The vanishing viscosity solution [6] of a Riemann problem for (3.1) is obtained by constructing a 
Lipschitz continuous map 

,R 



{su....s,)^Tl[T!'-\[...[Tl[u-]]]\ 



which is one to one from a neighborhood of the origin onto a neighborhood of v?" . Then we can uniquely 

^ and the wave sizes si, S2, ■ ■ ■ , sn such that 



determine intermediate states = ujo,oji, ■ ■ ■ ,ujn = u 
' is sufficiently small. 



1, 



u 




provided that \u^ 

By Theorem 3.2, each Riemann problem with initial date 

(3.7) uo 

admits a vanishing viscosity solution Uk, containing a sequence of rarefactions, shocks and discontinu- 
ities of the k-th family: we call Uk the fc-th elementary composite wave. Therefore, under the strict 
hyperbolicity assumption, the general solution of the Riemann problem with the initial data (3.5) is 
obtained by piecing together the vanishing viscosity solutions of the elementary Riemann problems given 
by (3.1)-(3.7). 

Indeed, from the uniform hyperbolicity assumption (3.2), the speed of each elementary k-th wave in 
the solution Uk is inside the interval [A/c_i,Aa;] if s <C 1, so that the solution of the general Riemann 



6 



STEFANO BIANCHINI AND LEI YU 



lairn solution 



erage matrix 



e: tilde all 



d:l 



e:tilde_r 



e:tilde_l 



problem (3.1)-(3.5) is then given by 



x/t < Ao 



(3.8) 



u{t,x) = I uk{t,x) Afc^i < x/t < Afc,fc = l,...,iV), 



x/t > A 



N- 



Remark 3.4. If the characteristic fields arc cither genuinely nonlinear or linearly degenerate, the admissible 
solution of Riemann problem (3.1)- (3.5) consists of N family of waves. Each family contains either only 
one shock, one rarefaction wave or one contact discontinuity. However, the general solution of a Riemann 
problem provided above may contain a countable number of rarefaction waves, shock waves and contact 
discontinuities. 

3.3. Cantor part of the derivative of characteristic for i-th waves. Recalling the solution (3.8) 
to the Riemann problem (3.1)-(3.5), we denote Xi{u^,u^) as the i-th eigenvalue of the average matrix 



(3.9) 



A(m^, w") = / AiOu'^ + (1 - 9)u'^)d0, 







and li{u ,u ), fi{u ,u ) are the corresponding left and right eigenvector satisfying li ■ fi = 5ij and 
\fj\ = 1, for every i,j G {!,..., N}. Define thus 

(3.10a) Xi(t, x) = Xi{u{t, X—), u{t, x+)), 



(3.10b) 



h{t,x) = r^{u{t,x-),u{t,x+)), 



(3.10c) li{t, x) = li{u{t, x—),u{t, x+)). 

Since the fi, U have directions close to r^, one can decompose D^u into the sum of N measures: 



N 



k=l 



where Vi = h ■ D^u is a scalar valued measure which we call as i-th wave measure [8]. 

In the same way we can decompose the a.c. part D'^u, the Cantor part D'^u and the jump part 
Di'^'^Pw of D^u as 



N N N 

fc=l k=l k=l 

We call the Cantor part of Vi and denote by 

the continuous part of Vi. According to Volpert's Chain Rule 

(3.11) D,X,{u) = \7X,{u){Dl'u + Dlu) + [A,(u+) - A,(w-)]4, 
and then 

(3.12) DlX,{u) = VX, -Dlu^Y. (^^' ■ ^^^H- 



.comp_laiiibda (3.13) 



We define the i-th component of DxXi{u) as 



[D.X,{u)], - (VA, • h)vr'' + [X^{u+) - K{u-)] 



and the Cantor part of i-th component of DxXi{u) to be 
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4. Main SBV regularity argument 

Following [7], the key idea to obtain SBV-like regularity for Vi is to prove a decay estimate for the 
continuous part of Vi. We state here the main estimate of our paper. 

Theorem 4.1. Consider the general strictly hyperbolic system (3.1), and suppose that the i-th character- 
istic field is genuinely nonlinear. Then there exists a finite, non-negative Radon measure nY^'^ on M.'^ x R 
such that for t > t > 

(4.1) I 

for all Borel subset B of 



vrHt)\{B) < o{i) 



C{B) 



+ ^if^\[t-T,t + T\ X K) 



Different from [7], we assume only one characteristic field to be genuinely nonlinear and no other 
requirement on the other characteristic fields. 

Once Theorem 4.1 is proved, then the SBV argument develops as follows [7]. 

Suppose at time t = s, Vi{s) has a Cantor part. Then there exists a £^-ncgligible Borel set K with 



(s)(/f) > and D'''"'Pv^{K) 



0. Then for all s > r > 0, 



0<\vM\iK) = \vr'is)\{K)<O{l) 



T 
,ICJ/- 



+ ^iI'^-'{[s-t,s + t] X M) 



Since C^{K) ~ 0, we can let t — > 0, and deduce that t^l'"''{{s} x M) > 0. This shows that the Cantor 
part appears at most countably many times because fiY-"^ is finite. 

Then, we can have the following result which generalizes Corollary 3.2 in [7] to the case when only one 
characteristic field is genuinely nonlinear and no assumptions on the others. 

Corollary 4.2. Let u be a vanishing viscosity solution of the Cauchy problem for the strictly hyperbolic 
system (3.1), and assume that the i-th characteristic field is genuinely nonlinear. Then Vi(t) has no 
Cantor part out of a countable set of times. 

As we see in the scalar case, by proving the SBV regularity of the solution under the genuinely nonlinear 
assumption of one characteristic field, we can deduce a kind of SBV regularity of the characteristic speed 
for general systems. 

Unlike the scalar case, we do not have the maximum principle to guarantee the small variation of u in 
the triangle T{t(),xo) defined in (2.2). However, in the system case, we have the following estimates for 
the vanishing viscosity solutions. 

For a < b and t > 0, we denote by Tot.Var.{u(r); ]a, 6[} the total variation of u{t) over the open 
interval ]a,b[. Moreover, consider the triangle 



I (i, x) : T <t < (b-a) /2t], a + 77^ < .t < - Tyt | . 



The oscillation of u over A^']J will be denoted by 

Osc.{u; Al-;}:=sup{\u{t,x)-u{t',x')\: {t,x), (t', x') G A^;','} . 

We have the following results. 

Theorem 4.3 (Tame Oscillation). [6] There exists C" > and fj > such that for every a < b and 
r > 0, one has 

Osc.{u; ^ C" • Tot.Var.{u(T); ]a,6[}. 

Adapting the proof of the scalar case, we can prove the main Theorem 1.1 of this paper: the proof of 
this theorem will be done in Section 7. 

5. Review of wave-front tracking approximation for general system 

To prove Theorem 4.1, we use the front tracking approximation in [3] which extends the one in [8] 
to the general systems. Since the construction is now standard, we only give a short overview about 
existence, compactness and convergence of the approximation, pointing to the properties needed in our 
argument: more precisely, we will only consider how one construct the approximate wave pattern of the 
A:-th genuinely nonlinear family (Section 5.1.2). 
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The main point is that, for general systems, the accurate/simphfied/crude Riemann solvers for the k-th 
wave coincides with the approximate/simplified/crude Riemann solvers when all families are genuinely 
nonlinear (see below for the definition of accurate/simplified/crude Riemann solvers). This means that 
the wave pattern pf the fe-th genuinely nonlinear family will have the same structure as if all other families 
arc genuinely nonlinear: by this, we mean that shock-shock interaction generates shocks, the jump in 
characteristic speed across fc-th waves is proportional to their size, and one can thus use the fc-componcnt 
of the derivative of Afc (3.13) to measure the total variation of Vk- 

5.1. Description of front tracking approximation. Front tracking approximation is an algorithm 
which produces piecewise constant approximate solutions to the Cauchy problem (3.1). Roughly speaking, 
we first choose a piecewise constant function Uq which is a good approximation to initial data uq such 
that 

(5.1) Tot.Var.{w^} < Tot.Var.{uo}, H^o ^ "oIIli < e, 

and Uq only has finite jumps. Let xi < ■ • • < Xm be the jump points of Uq. For each a = 1, . . . , m, we 
approximately solve the Riemann problem (see section 3.2, just shifting the center from (0,0) to {0,Xa)) 
with the initial data given by the jump [uff{xa—), [uQ(xa+)] by a function w(t,x) = (t>{jzj^) where is 
a piecewise constant function. The straight lines where the discontinuities locate arc called wave- fronts 
(or just fronts for short). The wave-fronts can prolong until they interact with other fronts, then at the 
interaction point, the corresponding Riemann problem is approximately solved and several new fronts 
are generated forward. Then one tracks the wave- fronts until they interact with other wave- fronts, etc... 
In order to avoid the algorithm to produce infinite many wave-fronts in finite time, different kinds of 
approximate Riemann solvers should be introduced. 

5.1.1. Approximate Riemann solver. There are two kinds of approximate Riemann solvers defined for 
interactions between two physical wave-fronts. Suppose at the point {ti,xi), a wave-front of size s' 
belonging to fc'-th family interacts from the left with a wave-front of size s" belonging to fc"-th family 
for some k' , k" G {1, ■ • • , N} such that 

Assuming that \u^ — vP-] sufficiently small. At the interaction point, the Riemann problem with the initial 
data data [u^ ^ vP] will be solved by approximate Riemann solver. 

• Accurate Riemann Solver replaces each elementary composite wave of the exact Riemann solution 
(refers to Uk in (3.8)) by an approximate elementary wave which is a finite collection of jumps 
traveling with a speed given by the average speed Afc given by(3.10a), and the wave opening (i.e. 
the difFcrcncc in speeds between any two consecutive fronts) is less than some small parameter e 
controlling the accuracy of the approximation. 

• Simplified Riemann Solver only generates approximate elementary waves belong to fc'-th and 
fc"-th families with corresponding size s' and s" as the incoming ones if k' ^ fc", and approximate 
elementary waves of size s' -\- s" belong to fc'-th family if k' = k" . The simplified Riemann solver 
collects the remaining new waves into a single nonphysical front, traveling with a constant speed 
A, strictly larger than all characteristic speed A. Therefore, usually the simplified Riemann solver 
generate less outgoing fronds after interaction than the accurate Riemann solver. 

Since the simplified Riemann solver produces nonphysical wave-fronts and they can not interact with 
each other, one only needs a approximate Riemann solver defined for the interaction between, for example, 
a physical front of the fc-th family with size s, connecting u'^ , u^ and a nonphysical front (coming from 
the left) connecting the left value u^ and u^^ traveling with speed A. 

• Crude Riemann Solver generates a fc-th front connecting u^ and u*^ = Tg'^[w^] traveling with 
speed Xi and a nonphysical wave-front joining u^'^ and u^, traveling with speed A. In the following, 
for simplicity, we just say that the non-physical fronts belong to the {N + l)-th characteristic 
field. 

Remark 5.1. We can assume that at each time t > 0, at most one interaction takes place, involving 
exactly two incoming fronts, because we can slightly change the speed of one of the incoming fronts if 
more than two fronts meet at the same point. It is sufficient to require that the error vanishes when 
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e ^ 0. 

To simplify the analysis, we assume that the fronts satisfy the Rankine-Hugoniot conditions exactly. 

5.1.2. The approximate Riemann solvers for genuinely nonlinear waves. If the fc-th characteristic family 
is genuinely nonlinear, the elementary wave Uk is either a shock wave or a rarefaction wave. The key 
example of the accurate Riemann solver is thus to consider how these two solutions are approximated. 
If fc-th elementary wave Uk in (3.8) is just a single shock, for example 

Uk = 



,R 



x/t < a, 
x/t > a, 



where a is the speed of shock wave, then the approximated fc-th wave coincides the exact one (apart from 
the speed in case, see the above remark). 

If Uk is a rarefaction wave of the fc-th family connecting the left value and the right value , for 
example, if := Tj^'[u'"] and 

{u^ x/t<\i{u^), 
T^.[u^ x/t e [K{u^),\,{v^% x/t = K{T^.[u^), 
x/t> \i{u^), 

where s* G [0,s]. Then the approximation Uk is a rarefaction fan containing several rarefaction fronts. 
More precisely, we can choose real numbers = sq < si < • • • < s„ = s. and define the points 
Wi := Tg. [v}"\^ i = 0, . . . , ri, with the following properties. 



w, 



i+l 



and the wave opening of consecutive wave-fronts are sufficiently small, i.e. 

a';[u^]{s,+i) ~ a^^[u^]{s,) < e, = 0, . . . , n - 1. 

where the function cr^' is defined in Theorem 3.2. We let the jump [oj^, oj^-i-i] travel with the speed 
di :— Afc(a;i, oji-i-i) (3.10a), so that the rarefaction fan becomes 

v}" x/t < (Tl, 

Uk ^ ■{uji a-i < x/t < CTi+i, I 1, . . . , n - 1, 

^ x/t>an. 

5.1.3. Interaction potential and BV estimates. Suppose two wave-fronts with size s' and s" interact. In 
order to get the estimate on the difference between the size of the incoming waves and the size of the 
outgoing waves produced by the interaction, we need to define the amount of interaction I(s', s") between 
s' and s". 

When s' and s" belong to different characteristic families, set 

I(s',s") = |s's"|. 

If s', s" belong to the same characteristic family, the definition of I(s',s") is more complicated (see 
Definition 3 in [3]). Wc just mention that if s', s" are the sizes of two shocks which have the same sign, 
traveling with the speed a' and a" respectively, then the Amount of Interaction takes the form 

(5.2) I{s',s") = \s's"\\a' -a"\, 

i.e. the product of the size of the waves times the difference of their speeds (of the order of the angle 
between the two shocks). 

To control the Amount of Interaction, the following potential is introduced. 

At each time t > when no interaction occurs, and u(t, ■) has jumps at xi, . . . , Xm, we denote by 



Wl, . . . , UJm, Si, . 



«1, 



5 ; 



their left states, signed sizes and characteristic families, respectively: the sign of Sq, is given by the 
respective orientation of dT!:[u]/ds and rk, if the jump at Xa belongs to the fc-th family. The Total 
Variation of u will be computed as 

V{t)^V{u{t)) ■.= Y.\s^\. 
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Following [5], we define the Glimm Wave Interaction Potential as follows: 



d:gp (5.3) 



Q{t) ^ QHt)) := 



E 



ia,=ifl<N+l 



Denoting the time jumps of the Total Variation and the Glimm Potential as 
AV{t) = V{t+) - V{t-), AQ(t) = Q(t+) - Q(t-), 
the fundamental estimates are the following (Lemma 5 in [3]): in fact, when two wave- fronts with size 



s', s" interact, 



(5.4a) 



AQ(r) = -0(l)Z(s',s"), 



(5.4b) AV{t) =0{l)I{s',s"). 

Thus one defines the Glimm Functional 

(5.5) T{t):^V{t) + CoQ{t) 

with Co suitable constant, so that T decreases at any interaction. Using this functional, one can prove 
that e- approximate solutions exist and their total variations are uniformly bounded (see section 6.1 of 
[3]). 

5.1.4. Gonstruction of the approximate solutions and their convergence to exact solution. The construc- 
tion starts at initial time t ~ with a given e > 0, by taking Mq as a suitable piecewise constant 
approximation of initial data uq, satisfying (5.1). At the jump points of Uq, we locally solve the Riemann 
problem by accurate Riemann solver. The approximate solution u'^ then can be prolonged until a first 
time ti when two wave-fronts interact. Again we solve the Riemann problem at the interaction point by 
an approximate Riemann solver. Whenever the amount of interaction (see Section 5.1.3 for the defini- 
tion) of the incoming waves is larger than some threshold parameter p ~ p{e) > 0, we shall adopt the 
accurate Riemann solver. Instead, in the case where the amount of interaction of the incoming waves 
is less than p, we shall adopt two different types of simplified Riemann solvers. And we will apply the 
crude Riemann solver if one of the incoming wave-front is non-physical front. One can show that the 
number of wave-fronts in approximate solution constructed by such algorithm remains finite for all times 
(see Section 6.2 in [3]). 

We call such approximate solutions e-approximate front tracking solutions. At each time t when there 
is no interaction, the restriction u'^{t) is a step function whose jumps are located along straight lines in 
the (t, a;)-plane. 

Let {ei/}^]^ be a sequence of positive real numbers converging to zero. Consider a corresponding 
sequence of e^-approximate front tracking solutions u'^ := u^" of (3.1): it is standard to show that 
the functions t i— > u^it^ •) are uniformly Lipschitz continuous in L}- norm, and the decay of the Glimm 
Functional yields that the solutions u'^{t,-) have uniformly bounded total variation. Then by Helly's 
theorem, u'^ converges up to a subsequence in h,l^^(R'^ x M) to some function u, which is a weak solution 
of (3.1). 

It can be shown that by the choice of the Riemann Solver in Theorem 3.2, the solution obtained by 
the front tracking approximation coincides with the unique vanishing viscosity solution [6] . Furthermore, 
there exists a closed domain D C L^(M, 51) and a unique distributional solution u, which is a Lipschitz 
semigroup 2? x [0,+cxd[— > T) and which for piecewise constant initial data coincides, for a small time, 
with the solution of the Cauchy problem obtained piecing together the standard entropy solutions of the 
Riemann problems. Moreover, it lives in the space of BV functions. 

For simplicity, the pointwise value of u is its representative such that the restriction map t u{t) 
is continuous form the right in and x i— S- u{t,x) is right continuous from the right. 

5.1.5. Further estimates. To each u", we define the measure p\, of interaction and the measure n]p of 
interaction and cancellation concentrated on the set of interaction points as follows. If two physical front 



fronts belonging to the families i,i' G {1, . . . , N} with size s' 



interact at point P, we denote 
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|s I + |s I - |s + s I 1 = 1, 


The wave size estimates (Lemma 1 in [3]) yields balance principles for the wave size of approximate 
solution. 

More precisely, given a polygonal region T with edges transversal to the waves it encounters. Denote by 
^i!tn' ^i^out the positive (+) or negative (— ) i-th waves in entering or exiting T, and let W^j^ = 



f.in 



■ w w . 

f ,111 ' f ,out 



* ' I/, out 



Qjjt . Then the measure of interaction and the measure of interaction- 



cancellation control the difference between the amount of exiting i-th waves and the amount of entering 
i-th waves w.r.t. the region as follows: 

(5.6a) IKout-l^:i„l<0(lK(r), 



e:bl_l 



de 


: of muic 




Ss_decay 



(5.6b) 



w. 



i± 



The above estimates are fairly easy consequence of the interaction estimates (5.4) and the definition of 



By taking a subsequence and using the weak compactness of bounded measures, there exit measures 
and ^i^'^ on R+ x M such that the following weak convergence holds: 

(5.7) ^il ^\ n 



ic 



IC 



5.2. Jump part of z-th waves. The derivative of is clearly concentrated on polygonal lines, being a 
piecewise constant function with discontinuities along lines. To select the fronts of converging to the 
jump part of u, we use the following definition. 

Definition 5.2 (Maximal (e'', e^)-shock front). [8] A maximal (e°, e^)-shock front for the i-th family of an 
Eiz-approximate front-tracking solution u'^ is any maximal (w.r.t. inclusion) polygonal line {t,^'^{t)) in 
the (t,x)-plane, to < t < ti, satisfying: 

(i) the segments of 7"^ are i-shocks of u'^ with size \s^\ > e°, and at least once \s'^\ > e^; 

(ii) the nodes are interaction points of u''; 

(iii) it is on the left of any other polygonal line which it intersects and which have the above two 
properties. 

Let M^^Q be the number of maximal (e°, e^)-shock front for the i-th family. Denote 



'^(eO,ei),m ' [*(£o',€i ) ,m ' ^(co',ei),m 



m ~ 1, 



as the maximal (e*^, e^)-shock fronts for the i-th family in u'^ . Up to a subsequence, we can assume that 



M, 



is a constant independent of ly because the total variations of m'^ are bounded. 



Consider the collection of all maximal (e°, e^)-shocks for the i-th family and define 



U Graph(7(";^^^i)_^), 



and let {e^JjfegN, {cfelfeeN be two sequences satisfying < 2'^e" < \ 0. 

Up to a diagonal argument and by a suitable labeling of the curves, one can assume that for each fixed 
k, m the Lipschitz curves j'f'l, l^ converge uniformly to a Lipschitz curve 7? i> • Let 

:=UGraph(7|,o.,i),,J. 



d:lnu 



d: vnujump 



denote the collection of all these limiting curves in u. 
For fixed (6^,6^), we write for shortness 



(5.8) 

and define 
(5.9) 



l'^{t,x) ■.= Uu''{t,x^),u''{t,x+)) 



V 



yjump 

i,(eO,ei) 



1" ■ n'^< „ , 
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Following the same idea of the proof of Theorem 10.4 in [8], the next lemma holds if only the i-th 
characteristic field is genuinely nonlinear. 

Lemma 5.3. The jump part of Vi is concentrated on ^\ 

Moreover there exists a countable set Q C IR+ x M, such that for each point 

P=(r,0 = (r,7™(r)) 

where i-th shock curve 7*„ is approximated by the sequence of {e^ , e^)-shock fronts j'^^l ^ of the approx- 
imate solutions , the following holds 



(5.10a) 



lim limsup 



\ 



sup 



\u'^{t, x) — u 



= 0, 



e:rightlim I (5.10b) 



ge of vijump (5.11) 



lim limsup 



sup 



kt'^(t, x) — 



V 



(t,x)l=B(J',r) 

Moreover, we can choose a sequence {vk}^^i such that 



N 



.jump 



weak* — lim > v^'', ' i,^"?'? 



goflefteigen (5.12) 



The key argument of the proof is that we can use the tools of the proof of Theorem 10.4 in [8] because 
the wave structure of the i-th genuinely nonlinear family has the following properties: 

(1) the interaction among two shocks of the i-th family generates only one shock of the fc-th family, 

(2) the strength of i-th waves can be measured by the jump of the i-th characteristic speed A,;, 

(3) the speed of i-th waves is very close to the average of the jump of across the discontinuity. 

These properties are a direct consequence of the behavior of the approximate Ricmann solvers on the 
i-th waves, if the i-th family is genuinely nonlinear (Section 5.1.2). 

Sketch of the proof. Let Q be the set defined by all jump points of the initial datum, the atoms of /i^'-' 
(see (5.7)). 

For any point P G 5^' \ O, if (5.10a) or (5.10b) does not hold, then this means that the approximate 
solutions u'^ has some uniform oscillation: indeed, because of the L^-convergence in M^, one can find 
points {t,x) arbitrarily close to P such that u'''{t,x) u{t,x), and the fact that the limits are not 
means that there are other points {t',x') arbitrarily close to P such that \u'^{t',x') — u{t,x)\, at least for 
a subsequence of v. The analysis of the proof of Theorem 10.4 in [8] shows that in this case the Amount 
of Interaction and Cancellation is uniformly positive in every neighborhood of P, and thus P is an atom 
of /x^*^, contradicting to P ^ 0. 

For P ^ U 0, if vf^™^{P) > 0, i.e. P is a jump point of u, by the similar argument of Step 8 
in the proof of Theorem 10.4 in [8] this shows that the waves present in the approximate solutions are 
canceled, and thus ^^^(P) > 0. It is impossible since P ^ Q. This concludes that i;^"™'' is concentrated 
on because by (5.10) the jumps in the approximate solutions are vanishing in a neighborhood of 
every P ^ U 6. 

We are left with the proof of (5.11). At jump point (i,7^^o ^i) ^ \ ®' according to (5.10a), 

(5.10b), there exist a sequence (t'^^j'^^l ^i-, „j(0(^'')) such that 



lim 



(t" 



and its left and right values converges to the left and right values of the jump in (t, 



,e-'-),m 



it))- 



^^,u^), 

and since its eigenvalues arc uniformly separated the same continuity holds for its eigenvalues Xk{u^, ^^)i 



Since / S C^, by the definition (3.9) the matrix A{u^, u^) depends continuously on the value {i 
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t:bde 



vibde 



ss needec 



:mu i nu est 



left eigenvectors lkiu^,u^) and right eigenvectors ffc(w^,M^). Using the notation (3.10a) and (5.8), one 
obtains 

(5-13) f,(<,7(,o,,.),„(0) =li^l^r(^^7("j,,I).™(^^)), 

and similar limits holds for f,;, A^. 

Up to a subsequence {i^fc}, from the convergence of the graphs of ,!^,'^o'^i\ to =5^* and (5.10a), (5.10b), 
it is fairly easy to prove that 



DuL. g-i = lim Du^L. 



According to (5.9), (5.13) and (5.14), one concludes the weak convergence of ir*;''^"^^^ to wf™^ 



□ 



6. Proof of Theorem 4.1 



6.1. Decay estimate for positive waves. The Glimm Functional for BV functions to general systems 
has been obtained in [5], and when u is piecewise constant, it reduced to (5.3): and we will write it as Q 
also the formulation of the functional given in [5]. Moreover, for the same constant Co > of the Glimm 
Functional T{t) (5.5), the sum Tot.Var.(M) + CoQ(u) is lower semi-continuous w.r.t the norm (see 
Theorem 10.1 of [8]). 

For any Radon measure ^, we denote [^J"*" and [^]~ as the positive and negative part of fi according 
to Hahn- Jordan decomposition. The same proof of the decay of the Glimm Functional T(t) yields that 
for every finite union of the open intervals J = Ii U ■ ■ ■ U Im 



.1) 



(J) + CoQiu) < liminf {K]±(J) + CoQiu'')} , i = 1, . . . , n. 



I/— f oo 

as w"^ — > u in L^. 

In [8, 9] the authors prove a decay estimate for positive part of the i-th wave measure under the 
assumption that i-th characteristic field is genuinely nonlinear and the other characteristic fields are 
either genuinely nonlinear or linearly degenerate. By inspection, one can verify that the proof also works 
(with a little modification) under no assumptions on the nonlinearity on the other characteristic fields, 
since the essential requirements of strict hyperbolicity and of the controllability of interaction amounts 
by Glimm Potential still hold: the main variation is that one should replace the original Glimm Potential 
in [8] with the generalized one given in [5]. 

We thus state the following theorem, which is the analog of Theorem 10.3 in [8]. 

Theorem 6.1. Let the system (1.1) be strictly hyperbolic and the i-th characteristic field be genuinely 
non-linear. Then there exists a constant C" such that, for every < s < t and every solution u with small 
total variation obtained as the limit of wave- front tracking approximation, the measure [ui(t)]'^ satisfies 

-iB) + [Qis)-Qit)] 



(6.2) 

for every B Borcl set in 



[iH{t)V{B) < C" 



t 



The estimate (6.2) given half of the bound (4.1). 

6.2. Decay estimate for negative waves. To simplify the notation, we omit the index (6*^,6^) in 
^r'(eO™i) ^^^^ proof. In order to get the uniform estimate for the continuous part 11^"''^°"' := 

vf — ■;;'^'J"™P -we need to consider the distributions 



Mr :=9t< + 5.(AX), 



i^,jump 



6.2.1. Estimate for /i.^. Let : [t^,t+] — >■ K, 
graphs are the i-th wave-fronts of m"^ and define 



be time-parameterized segments whose 



— u{t,y^{t)- 



u{t,yyn{t)+), te]T„ 



For any test function £ C^{1 



X 



on obtains 



(6.3) 



/ 



>dfj.\ 



+ xl 



E 

771 — 1 



[0(^m,2/m(T+)) - (j){T^,y,niT^)]h ' («„ " "„)• 
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and strength 



action of pw 



For any m, since the i-th characteristic field is genuinely nonlinear, one has 

\k{u'^,iL''')-k{u'^)\=Oil)\ul-ut\, 
where u^^ = T"]. [u^J for some size Si. Then it follows from (3.6) that 

(6.4) - • (uf„ - ut). 

Let {{tk,Xk)}k be the collection of points where the i-th fronts interact. The computation (6.3) yields 
that yuj' concentrates on the interaction points, i.e. 

k 

where pk is the difference between the strength of the i-th waves leaving at {tk,Xk) and the i-th waves 
arriving at (tk ,Xk)- We estimate the quantity pk depending on the type of interaction: 

Since in [7], it is proved that the total size of nonphysical wave- fronts arc of the same order of e^, 
when decomposing u!^, we only consider the physical fronts. If at {tk-,Xk), two physical fronts with i-th 
component size sj;, s" interact and generate an i-th wave or a rarefaction fan with total size Si = s™, 
from (6.3) and (6.4), one has 

(6.5) Pk ^ -s[~ s". 

Notice that s' or s" or both may vanish in (6.5) if one of incoming physical fronts does not belong to the 
i-th family. 

According to the estimate in [3] (Lemma 1), the difference of sizes between the incoming and outgoing 
waves of the same family is controlled by the Amount of Interaction (see Section 5.1.3), so that one 
concludes 

\titmtk.xk)})<o{i)i{s.,,s';) 

and thus 

\K\{{tk}xm<0{l){T^{t-,)-T^{tt)). 

This yields 

(6.6) \^i';\{m^ y.^)<0{\)T''{Q), 
i.e. l/ij'l is a finite Radon measure. 



6.2.2. Estimate for /.i 



Let 7* 



771=1, 



^lf_o el) J be the curves whose graphs are 



the segments supporting the fronts of u'^ belonging to =3^^o and write 



For any test function G C^{^+ 

(6.7) - / 

which yields 



Ml 



by direct computation one has as in (6.3) 



E 

m— 1 



e:q_k 



where (t^, Xk) are the nodes of the jumps in ^^^o^^ij and the quantities qk can be computed as follows: if 
the i-th incoming waves have sizes s' and s", and the outgoing i-th shock has size s, then (see [7]) 



'.8) Qk 



{tk,Xk) terminal point of a front not merging into another front, 
(^fe, Xk) initial point of a maximal front, 
- s" (tk, Xk) merging point of two fronts, 

[tk, Xk) interaction point of a front with waves not belonging to ^(^o ^ly 

Except in the case when (r^, Xk) is the terminal points of the front of =^^o'gi) which ends without merging 
into another, one has by the interaction estimates 

Qk < fl]P{Tk,Xk). 
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balaiice_jump 



balance 



In fact, since s < on shocks the second case of (6.8) imphes qk < 0. 

Suppose now that (r^, Xk) is a terminal point of an (e", e^)-shock front 7„i. By the definition of (e°, e^)- 
shock, for some t < Tk the shock front 7™ has size sq < — e^, and at {Tk,Xk) the size si of the outgoing 
i-th front must be not less than —e^ as a result of interaction-cancellation among waves. Hence we obtain 



ki|< 0(1^^(7.0 



This yields 



qk ^ ~si + (si + qk) < ^— ^(e' - e") + Oil)^ll{tk,Xk) < ^^/ii^(7fc) + Oil)^ilitk,Xk). 



Since the end points correspond to disjoint maximal z-th fronts, due to genuinely nonlinearity, it follows 
that 

qk< 0(l)Mi''(M+ X M), 

{ik-,^k) end point 

SO that it is a uniformly bounded measure. We thus conclude that the distribution 

{tk^Xf^) end point 

is non-negative, so it is a Radon measure and thus also ^^'■'""^p is a Radon measure. 
In order to obtain a lower bound, one considers the Lipschitz continuous test function 



</'a(i) := X[0,T+a](0 



t-T 



X[T,T+a]{t), a > 0, 



which is allowed because is a bounded measure. Since is non-negative, one obtains 



/2'^([0,T] xR) < 



ic 



qk4>a{tk) 

{tk TXk) end point 



< 



< 



+ xl 



T+a 



O(iy,C([0,r + a] xM). 



Letting a \ and since [Uj'''^"™P(M)](0) is negative, one concludes 

/2"([0,r] xE) < -[i;,"'j"'"P(r)](E) + O(l)/iiC([0,T + a] x K) < O(l)T^(0). 

We conclude this section by writing the uniform estimate 

-O(l)T"(0) < < 0{l)fji^. 

In particular, the definitions of the measures /i^, ^^'J"™p give the following balances for the i-th waves 
across the horizontal lines: 



(6.9a) 



[vUt+)] (M) - [<(t-)] (K) = ({0 X M), 



balance jump (6.9b) 



[<'j""^p(t+)](M) - K'j""^p(t-)](E) = Mr'""^({n X 



The limits are taken in the weak topology. Notice that we can always take that t n- v'^Uiit), v'('^^'^^{t) is 
right continuous in the weak topology. 
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6.2.3. Balances of i-th waves in the region bounded by generalized characteristics. We recall that a min- 
imal generalized i-th characteristic is an absolutely continuous curve starting from (to,a;o) satisfying the 
differential inclusion 

x''it;to,xa) ■.^iJim{x''{t):x''{to)=xa, (t) £ [A,(w"(t, .t(0 + ) , A, a;(i)-))] } 

for a.e. t > to. Given an interval / = [a, 6], we define the region ^['^J^^j''^^ bounded by the minimal i-th 
characteristics a{t), b(t) of starting at {to, a) and {to,b) by 



A 



a,b] 



\{t,x) : to < t < to + T, a{t) <x< 6(t)|, 



and its time-section by I{t) 
intervals {Ii}fLi, and set 



i{t), b{t)]. Let J := /i U /2 U • • ■ U Im be the union of the disjoint closed 



M 



J(t) :=/i(t)U---U/M(i), 



U4 



flux control 



e:vibalance 



ijumpbalance 



We will now obtain wave balances in regions of the form A'^j^*°'^\ Due to the genuinely non-linearity of 
the z-th family, the corresponding proof in [7] works, we will repeat it for completeness. 

The balance on the region has to take into account also the contribution of the flux ^J' across 

boundaries of the segments due to the definition of generalized characteristic and the wave- front 

approximation, it follows that is an atomic measure on the characteristics forming the border of 
and moreover a positive wave may enter the domain A'^j'^*°''^'' only if an interaction occurs at 
the boundary point {t, x), which gives the estimate 

(6.10) $r({(t»}) <0(lVf 
One thus obtains that 

(6.11) [<(t)] (J(r)) - [<(to)] (J) = t^-iA"/""^^) + ^•^{A''/'"-^) + 0{l)e,, 

where the last term depends on the errors due to the wave-front approximation (a single rarefaction front 
may exit the interval Im at to). 

The same computation can be done for the jump part ^j^^-J"™?^ obtaining 



.12) 



[<^J"'"P(J(t))] (r) - [v'""''(io)] (J) = ;,^'J"-P (a;;' (*"•"') + $J'd"-P(A^'l'0'-)) 



.l',{to,T}\ 



Since the flux <j)^'J"™P only involves the contribution of (e", e^)-shocks, it is clearly non-positive. 
Subtracting (6.12) to (6.11), one finds the following equation for ■yj^''^"'^*: 



K'™"*(r)] ( j(r)) - K^™"*(to)] (J) = (A^r - 1^) (^r) + (*r - r^™p) (a 

Denote the difference between the two fluxes by 



J 



0(l)e.. 



Since (i)^'J"™P removes only some terms in the negative part of one concludes that 

for phicont I (6.13) $r - < [$,r]^ < m!.^- 

Setting 

and using the estimate l^fl 5; C(1)mI?: has 
for muicont | (6.14) - /<'j™P < 0(1)^^,^. 
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6.2.4. Decay estimate. Due to the semigroup property of solutions, it is sufficient to prove the estimate 
for the measure [w^"''^""' (t = 0)]^. Consider thus a closed interval / = [a,b] and let z{t) := a{t) — b{t) 
where 

a{t) := x''{t;0,a), b{t) := {t;0,b) 
and the minimal forward characteristics stating at t = from a and b. For £^-a.e. t one has 

z{t) = X,it,b{t)) ~ ~X,it,a{t)). 

By introducing a piecewisc Lipschitz continuous non-decreasing potential $ to control the waves on the 
other families [8], with $(0) = 1, one obtain 



(6.15) 
where 



at) ■= {kit, ait)) - Kit, ait)-)) + (A,(i, bit)+) ~ Kit, bit))) . 

We consider two cases. 
Case 1. If 



iW-$W^W< JK(o)](/) 



for alH > 0, then 
d 



dt 



- /' <E>(r)£ir 

"/o *W'i-,(t) ^ e-/o - mzit)] <"—- [<(0)](/) < O(l)[<(0)](/). 



Integrating the above inequality from to r and remembering that $(0) — 1, one has 

-C'il) = z(0) < e-J^'-i *(")'^"z(t) - z(0) < 0(1)t[<(/)] < 0(1)t[ 
which reads as 

CHI) 



«r™"*(o)](/), 



\ii)<oii) 



Case 2. Assume instead that 



zW-$(t)z(t)>-[<(0)](/) 



at some time t > 0. By the inequality (6.15) and the balance (6.11) one obtains 



i[<(0)](/)<i(t)-<i>(t)z(i) 



(6.16) 

Hence 
(6.17) 



< [<(t)](/(t))-$(t) + 0(l)e. 

< [<(o)] (/) + + *r (^FaS") - ^(^) + ^(i)^- 



K(o)](/)> 



i',(0,«) 



From (5.9) and the fact that the fronts in J^^i, satisfy Rankine-Hugoniot conditions (up to a negligible 
error), we have 



\t, ait)) = \,it, ait)+) - \,it, ait)-). 



and since 



\iu^) + Kiu^) 



< ^^(i)^^-^^!'. 



by the balance (6.12), we conclude that 



m > l[[v-:"it)] iait)) + [«r^"'"P(i)] ibit)) 2ei 



.18) 



> - 

- 4 

3 

> - 

- 4 



K'J"-P(t)](/(t))-2ei 

[<j-p(o)] (/) + A^r^"-p(A-;j*)) + *)) - 2.^ 
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Substituting (6.18) into (6.17), using the estimates (6.13), (6.14) and integrating (6.16), we obtain 

- °"(0)] (/) < 0(1)|^ + m!.^J(^4^) +e'+ 

This gives the estimate (4.1) for the case of a single interval for the approximate solution. 

By repeating the analysis in the case of a finite union of intervals, one obtains the same bound as 
above, and since is a Radon measure, the same result holds for any Borel sets, i.e. 

- [.r"*(0)] {B) < 0(l)|^ + ^i^-\A^) +e'+ ..|, 
where B is any Borel in M and 

:= {(T,a;''(T;0,a;o)) : 2; £ B, < r < 

As the solution is independent on the choice of the approximation, we can consider a particular 
converging sequence {u'^}u>i of Ciy-approximate solutions with the following additional properties: 

(6.19) Q{u''{0,-))-^Q{uo). 

By lower semi-continuity of [t;i(0)]~ + CoQ(u(0)) (6.1), one gets 

(6.20) K(0)]- +CoQ(ii(0)) < weak* -liminf {[<(0)]~ +CoQ(w"(0))} . 

From (6.20) and (5.11), up to a subsequence, one obtains for any open set _B C M, 

< hminf {[<(0)] "([/)+ C7oQ(u'^(0))} -CoQ(u(0))+ lim K'j™P(0)] (f/) 
= liminf{[<'™"*(0)] "([/) + CoQ(u'^(0))} -CoQ(u(0)) 

(6.21) 

< ]irnmiO{l)l^^+t,f'''(A^)+e'+e, + Qiu-^m ~ Q(u(0))l 

<O(l)|^+MP([0,i]xR)|, 

where //J^'' is defined as weak*-limit of measure /i,^'^^'' (up to a subsequence). Then the outer regularity 
of Radon measure yields the inequality for any Borel set. 
The above estimate together with Theorem 6.1 gives (4.1). 

7. SBV REGULARITY FOR THE i-TH COMPONENT OF THE i-TH EIGENVALUE 

This last section concerns the proof of Theorem 1.1. adapting the strategy of Section 2. 

Proof of Theorem 1.1. As in the scalar case, we define the sets 

Jr := {x e M : u^{t,x) ^ ii"(T,.T)}, 

Fr := {xeR: VX^{u{T,x)) ■ r,{u{T,x)) = O}, 
C := {(r,0 e ]R+ X E : ^ e Jr U Fr} , a:=J^UJV- 
By definition of continuous part 

\vr\Jr)\{r) = 0, 

and since 

VA,(u(t,F, \ J,)) ■r,(u(r,F,\ J,)) =0, 

we conclude that 

|VA,(u) • r,(«)<°"*(r)|(a) = \VKiu) ■ r,(u)«r"*(T)| ( J,) + |VA,(u) • r,(^/)<°"'(T)| (F, \ ,/,) = 0. 
For any (to^xo) e M+ x M \ C, there exist strictly positive bo ~ 5o(a;o,to)j cq = co(xo,to) such that 

|VAi • ri{u{to, x))\ > Co > 
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e :uni_conv 



e : local 



ex : nonsbv 



for every X in the open interval /o :—]— bo + xo,xo + bo[, because u(t(),x) is continuous at ^ Ctg- Hence 
by Theorem 4.3, we know the there is a triangle 

To := : \x - xo\ < b'^ - ri{t - to), < < - to < b'^/f]^ 

with the basis Iq :—] — b^ + xo,xo + 69 [c /q, such that 



(7.1) 



\VX^■r,{u{to,x))\ > y > 0, 



by taking b'^ ^ 1 in order to have that the total variation remains sufficiently small. 
Since wltq coincides with the solution to 

dtw + f{w)x = 0, 

Uto{x) \x - Xol < &0: 



(7.2) 



w{x,to) 



. Ixo-b'° "to {y)dy \x-xo\>bo, 

and by taking 69 sufficiently small, we still have that (7.1) holds for the range of w. In particular w is 
SBV outside a countable number of times, and the same happens for u in To- 

As in the scalar case, one thus verifies that there is a countable family of triangles {Ti}°^^ covering 
the complement of C outside a set whose projection on the t-axis is countable. The same computation 
of the scalar case concludes the proof: for any r chosen as in (2.3) 

\D''K-n\{R) < \D''X,-r,\{Cr) + \D''X, -n\(\jT,n{t = T}] =0. 

□ 

Similar to the scalar case, it is easy to get the following corollary from the Theorem 1.1. 

Corollary 7.1. Let u be the vanishing viscosity solution of the problem (1.1), (1-2). Then the measure 
{DuXi{u) ■ ri(u))(li{u) ■ Ux) has no Cantor part in R"'' x M. 

Remark 7.2. Consider the following equations 

Uh=0, 

\vt + iil + v + u)v)x =0. 

Since DxX2{{u^ v)) = DxU + 2DxV, then it is clear that 0^X2 can have a Cantor part, while from Theorem 
4.1 the second component 

2 

(DuX2 ■ 7-2) (I2 ■ {u, v)x) = — -r- (vux + {1 + U + 2v)vx) 

1 + u + Zv 

has not Cantor part. 
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